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^^0 ■ We discuss the problems involved in extracting the nuclear equation-of-state from 

' heavy-ion collisions. We demonstrate that the equation of state becomes effectively 

1-^' softer in non-equilibrium and this effect is observable in terms of collective flow 

^ ' effects. Thus, non-equilibrium effects must be included in transport descriptions 

^•fj , on the level of the effective mean fields. A comparison with transverse momentum, 

rapidity, and centrality selected flow data show the reliability and limitations of 
the underlying interaction which was derived from microscopic Dirac-Brueckncr 
(DB) results. 

1—i ■ 

Jm^ • Heavy-ion collisions open the possibility to explore the nuclear equation- 

p^ I of-state (EOS) far away from saturation. In particular, the collective flow is 

f^ ' connected to the dynamics during the high density phase of such reactions and 

0^ I thus yields information on the nuclear EOS EJ. The different components of 
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>!< • collective flow, m particuJaii, the in-plane and out-of-plane flow □ has recently 

^^ I attracted great interest dBij, since also the energy, centrality and transverse 

momentum dependence has been measured. 

From the theoretical point of view, however, the determination of the 
EOS, which is an equilibrium concept, is not straightforward in heavy ion 
O I collisions, which are highly non-equilibrium pmcesses. An appropriate start- 

2 ' ing point are the Kardanoff-Baym equations u for the non-equilibrium real- 

time Green functions. To arrive at transport descriptions several approxi- 
mations are usually introduced: the semi-classical approximation neglecting 
non-locality and memory effects, the quasi-particle approximation neglecting 
«2j ' the finite width of the spectral functions, and the T-Matrix approximation, 

C^ I specifying the self-energies in the Brueckner ladder approximation. Severjal 

attemps have recently been made to improve the first two approximations o. 
They are expected to be important e.g. for subthreshold particle production, 
but it is not known whether they are also relevant for flow observables. 

With the above approximations one arrives at a BUU-type transport 
equation with DB self-energies which should be specified for the general non- 
equilibrium configuration. Since these cannot be calculated one further in- 
troduces a local nuclear matter approximation neglecting spatial variations 
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Figure 1. Left: Time evolution of the quadrupole moment of the energy-momentum tensor 
Qzz (bold lines) and of the central density (in units of psat , dashed lines) at the center in 
central Au + Au collisions at 0.6 and 6.0 A.GeV. Right: Energy per particle of a two- 
Fermi-sphere configuration (CNM) (minus the relative kinetc energy) as a function of total 
density. This "reduced" EOS is shown for DB self energies for several relative velocities 
and for two (equilibrium) Skyrme parametrizations. 



and adopting a model for the local nuclear matter phase space distribution. 
Two models have mainly been used: One is the Local Density Approximation 
(LDA) assuming equilibrated nuclear matter at the local total density, usu- 
ally supplemented with an empirical momentum dependence taken from the 
real part of the optical model. The other is the Local phase space Configu- 
ration Approximation (LCA) which assumes a configuration of two interpen- 
etrating streams of nuclear matter with a given relative velocity (also called 
colliding nuclear matter, CNM), i.e. a two- Fermi-sphere configuration with 
Pauli-effects taken into account. In the first case this is directly the nuclear 
EOS. For the second case the self energies have not been calculated in the DB 
approach (but for the non-relativistic caseQ). However, we have developed a 
self-consistent extrapolatiqa starting from nuclear matter DB calculations as 
discussed in detail in ref. U. This then gives the model non-equilibrium self 
energies as a function of the configuration parameters, namely the two Fermi 
momenta and the relative velocity. 

The relevance of such non-equilibrium effects in the dynamical situation of 
heavy ion collisions is demonstrated in Fig. |l] (left side). Here the quadrupole 

^_^ Q^7n33 rj^ 1 1 rTi'2'2 

moment of the energy-momentum tensor Qzz — x'i'i+T^^+T'^'^ ^^ shown to- 
gether with the baryon density as a function of time. It is clearly seen that 
the local phase space is significantly anisotropic during times which are compa- 
rable with the compression phase of the collision independent of beam energy. 
Hence, non-equilibrium effects are present during a large part of the compres- 
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sion phase where one wants to study the EOS at high densities. Therefore 
non-equihbrium effects should be considered at least on the level of the ef- 
fective in-medium interaction, i.e. the mean fields in transport calculations 
should be taken for non-equilibrium configurations instead for equilibrated nu- 
clear matter. Actual transport calculations have verified that the anisotropic 
phase space iri. heavy ion reactions can be parametrized reliably by CNM 
configurations □. 

In order to have a qualitative feeling for the non-equilibrium effects we 
define an effective, "reduced" EOS for the CNM configuration. This is given 
by subtracting from the binding energy per particle for the CNM configuration 
with a given relative velocity Vrei the irrelevant kinetic energy of relative 
motion between the two streams. This quantity is shown on the right side of 
Fig. |l| for symmetric CNM configurations {kpi = kp^ ) for relative velocities 
Vrei = 0.8(0.99) corresponding to beam energies of 0.6(6.0) GeV per nucleon 
as a function of the total density. Vrei — corresponds to the isotropic case 
(equilibrated nuclear matter, i.e. one Fermisphere). As a reference we also 
show the widely usejd Skyrme parameterisations with a soft/hard EOS with 
K = 200/380 MeV 0. 

It is seen that the effective EOS becomes more attractive, i.e. softer, 
compared to the equilibrium EOS and saturates at a higher total density. 
This effect is understood qualitativelyin the following way: If there was no 
interaction between the two nuclear matter streams then the reduced energy 
per particle of the two-Fermi-sphere configuration were just the energy per 
particle of one Fermi sphere at half the total denisity. The minimum of the 
reduced EOS would then be shifted to twice psat ■ The interaction energy per 
particle can be estimated as twice the real part of the optical potential at the 
relative momentum which can be estimated as twice the Fermi momentum 
at the corresponding density. For saturation density this is attractive but 
decreases with increasing momentum. Thus we expect a minimum of the 
reduced EOS between one and twice Psat and below the saturation energy of 
equilibrium nuclear matter. This is what is seen qualitatively in Fig. 1. 

The question of interest is now whether these non-equilibrium effects are 
observable in terms of collective flow. We have therefore performed transport 
calculations using CNM self enpgies derived from DB mean fields calculated 
by the Tiibingen-group (DBF) t2l in the two approximations: the LDA which 
refers to equilibrated nuclear matter (DBF/LDA) , and the local configuration 
approximation (LCA) where the non-equilibrium DBF mean fields are used 
in the transport calculations (DBF/LCA). 

We discuss azimuthal distributions as function of transverse momen- 
tum, rapidity, centrality and beam energy. Such distributions have been 
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Figure 2. Flow per particle Fy/A and elliptic flow V2 at mid-rapidity (lAy'"'] < 0.1) 
m semi-central Au + Au reactions from SIS to AGS energies. The transport calculations 
til with microscopic DB mean flcjlf^ in the LCA (solid) and LDA (dashed) approaches are 
compared to the data (symbols) | 



parametrized in terms of a Fourier series dN/d(j) = ao(l + 2aicos(0) -I- 
2a2Cos(2(j))) with (f) the azimuthal angle with respect to the reaction plane. 
The coefficients ai and a2 describe the in-plane and out-of-plane components 
of collective flow, respectively, and they depend on transverse momentum 
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(PM)/(FP™JMproj), rapidity y(o) = Ycm/YP^^"^ , centrahty and 



energy. The in-plane or directed flow vi can also be characterized by the 
quantity Fy 
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as the slope of the mean in-plane trans- 
y<")=o 
verse momentum at mid rapidity. V2 is often referred to as the elliptic flow. 

Its sign describes the transition from out-of-plane (w2 < 0) to in-plane flow 
{v2 > 0) Q. The squeeze-out ratio Rn is connected to the elliptic flow through 
Rn = (1 — 2f2)/(l + 2w2)- Fig- shows the observables Fy/A and V2 for 
Au + Au collisions from SIS to AGS energies. Rather complete data are 
available from the EOS and FOPI collaborations, which ape. however, not 
completely consistent for V2. It is seen that the calculations Eil reproduce the 
flow Fy only for SIS energies, while at AGS energies the calculations oversti- 
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Figure 3. Rapidity distributions for Au + Au collisions at 400 AMeV selected according 
to bins of ERAT which is a measure of centrality (increasing with increasing centrality). 
Calculations are shown for the LCA (solid) and H)^ (dashed) approximations for the self 
energies and compared with the data (symbols) ri . 



mate the data. This is to be expected because the DB self energies do not 
reproduce the empirical saturation of the real part of the optical potential at 
energies above about 1 GeV, but rather lead to linearly increasing repulsion. 
This situation also appears in non-relatipstic descriptions using momentum 
dependent Skyrme-paramctrizations (see E3). The difference between the LDA 
and LCA descriptions is not significant because the DBF self energies are not 
strongly momentum dependent at saturation density. 

The description of the elliptic flow is generally much better because as a 
ratio it is not so strongly dependent on the absolute transverse flow. We see 
that both models generally reproduce the qualitative trend of the data, i.e. 
the transition from small V2 at 0.2 A. GeV to a preferential out-of plane flow 
(v2 < 0) which is maximal at 0.4 A. GeV with the subsequent transition back to 
in-plane flow around 4 A. GeV. There is now, however, a pronounced difference 
between the two approximations. The calculation without non-equilibrium 
effects (LDA) shows a larger squeeze-out effect and also has a larger transition 
energy from a out-of-plane to an in-plane emission of participant matter. This 
effect is reduced for the LCA calculation which can be understood by the 
effective softening of the equation of state during the initial non-equilibrium 
phase of the collision as discussed in connection with Fig. H. We also mention 
that the observed difference is of the samp_magnitude as the difference between 
a soft and hard EOS (see Fig. g of Ref. E3). Furthermore the inclusion of the 
non-equilibrium effects considerably improves the agreement with the data 
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which means that the underlying effective interaction (DB) is able to describe 
the dynamics of heavy ion collisions. 

In Fig. H we show rapidity distributions for Au + Au collisions at 400 
AMeV selected according to a quantity ERAT, the ratio between transverse 
and longitiidinal total momentum, which has been shown to be a measure of 
centrality 113. The rapidity distributions are reasonably well reproduced, in 
particular for central collisions (large ERAT). A general tendency is that the 
spectator particles observed at F*^") ~ 1 in the experiment are more stopped 
in the calculations. This effect is larger in LDA than LCA, because the field 
is more repulsive in the former case (see Fig. 1). 

To summarise the effective equation-of-state probed by the compression 
phase in energetic heavy ion reactions is in a significant way governed by 
local non-equilibrium. The anisotropy of the phase space lowers the binding 
energy per particle and makes the effective EOS seen in heavy ion reactions 
softer. This fact is reflected in the behaviour of collective flow phenomena. We 
conclude that "geometric" phase space effects should be taken into account on 
the level of the effective interaction when conclusions on the equilibrium EOS 
are drawn from heavy ion collisions. The comparison with experiments shows 
that microscopic many-body methods methods that are successfuU for nuclear 
structure also explain heavy ion reactions quantitatively at SIS energies Ej. 
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